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A new approach toward geometrical concept of black hole thermodynamics 
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Motivated by the energy representation of Riemannian metric, in this paper we study different 
approaches toward the geometrical concept of black hole thermodynamics. We investigate thermo¬ 
dynamical Ricci scalar of Weinhold, Ruppeiner and Quevedo metrics and show that their number 
and location of divergences do not coincide with phase transition points arisen from heat capacity. 

Next, we introduce a new metric to solve these problems. We show that the denominator of the 
Ricci scalar of the new metric contains terms which coincide with different types of phase transi¬ 
tions. We elaborate the effectiveness of the new metric and shortcomings of the previous metrics 
with some examples. Furthermore, we find a characteristic behavior of the new thermodynamical 
Ricci scalar which enables one to distinguish two types of phase transitions. In addition, we gen¬ 
eralize the new metric for the cases of more than two extensive parameters and show that in these 
cases the divergencies of thermodynamical Ricci scalar coincide with phase transition points of the 
heat capacity. 


I. INTRODUCTION 

One of the succinct semi-classical approaches for investigating the quantum nature of gravity is via black hole 
thermodynamics in AdS spacetime which is dual to that of a field theory in one dimension fewer [I|. Besides, an 
interesting development in the black holes studies comes from the fact that black holes are akin to thermodynamical 
system which can be essentially described by the laws of thermodynamics re-expressed in terms of properties of black 
holes [l-Q. 

Scientific researches in the black hole physics have presented a very deep and fundamental relationship between 
quantum gravity, holography and thermodynamics. Among these researches, the special interest is attracted by AdS 
black holes because some of them lead to the very interesting phenomenon called Hawking-Page phase transition Q . 
Phase transition plays an important role in order to explore thermodynamical properties of a system near the critical 
point. Usually, phase transitions are denoted by a discontinuity of a state space variable, specially heat capacity Q. 

Motivated by large applications of the geometrical concept of thermodynamics in the black hole phase transition, 
in this paper, we investigate strengths and shortcomings of different approaches toward the matter. First attempt 
was done by Weinhold [9| which introduced a metric on the space of equilibrium states where its components are 
given as the Hessian of the internal energy. Then, Ruppeiner introduced a metric which is defined as the negative 
Hessian of entropy with respect to the internal energy and other extensive quantities of a thermodynamical system 
0- It was shown that these two metrics are conformally equivalent to each other where the temperature is the 
conformal factor [llj |. Recently, Quevedo [l^ proposed a Legendre invariant metric, in which solved some of problems 
in Weinhold/Ruppeiner methods. 

The basic motivation for considering the geometrothermodynamics comes from the fact that this formalism describes 
in an invariant way the thermodynamic properties of a given thermodynamical system in terms of geometric structures. 
Although extracting phase transition in terms of curvature singularity is the main reason to consider the geometrical 
approach in thermodynamics, there are several examples in which the curvature singularities of the known metrics 
(Weinhold and Ruppeiner metrics and their Legendre transformations) are not located at thephase transition points 
and they have number of singularities in which are before/after the phase transition points [l3l Il4|. 

In this paper, in order to solve the mentioned problems, we introduce a metric where its curvature singularities are, 
exactly, located at the phase transition points. As we see later, this effective metric has different structure from the 
other known metrics. For this claim, we consider the black hole solutions with Born-Infeld (BI) and Maxwell sources 
in three and four dimensions and investigate phase transitions in various geometrothermodynamical methods. 
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II. HEAT CAPACITY 


In order to investigate the local stability of a black hole, one can use various methods related to different ensembles. 
In principle thermal stability can be carried out in the grand canonical ensemble by finding the determinant of the 
Hessian matrix of M{Xi) with respect to its extensive variables Xi [Hig. For static charged black holes, we usually 
regard the mass M as a, function of the entropy S and the charge Q. The number of thermodynamic variables 
depends on the ensemble that is used. The most well known approach for studying phase transition is in the context 
of canonical ensemble, hence the heat capacity of systems. Depending on the number of extensive parameters, the 
behavior of Hessian matrix is highly sensitive |17| and therefore most of physicists prefer to work in canonical ensemble. 
In this ensemble, the positivity of the heat capacity is sufficient to ensure thermal stability. In addition the system is 
considered to be in fixed charged and the heat capacity has the following form 




Ms 

Mss' 


( 1 ) 


where Ms = {^§-)q and Mss = are regular functions. Now we regard two different types of phase 

transitions. In type one, the changes in signature of the heat capacity is representing a phase transition. In other 
words, the roots of the heat capacity in this case are representing phase transition points which means one should 
solve Ms = 0. Phase transition concerning to the divergency of the heat capacity is denoted by type two. It means 
the singular points of the heat capacity are representing the phase transitions. Thus, we should consider Mss = 0 to 
obtain the phase transition of type two. 

In order to have a fitting geometrical approach for studying phase transitions, the thermodynamical Ricci scalar 
(TRS) must diverges in the mentioned both types of phase transitions. In what follows, we present a brief review study 
of several geometrical approaches with their shortcomings, and then, we propose a new effective metric concerning 
this issue. 


III. WEINHOLD METRIC 

In Weinhold method, one is considering appropriate extensive parameters such as entropy, electric charge and 
angular momentum, and their related intensive quantities such as temperature, electric potential and angular velocity 
with the mass of the black holes as a potential. The Weinhold metric is given by @ 

= MgZdX<^dX\ ( 2 ) 

where = d'^M (A'^) jdX'^dX^ and also A“ = A“ (S', iV®), where iV® denotes other extensive variables of the system. 
Considering a static charged black hole, one can find the following expression for the denominator of Weinhold Ricci 
scalar 


denomiTZw) = {MssMqq - Ml^f M^ (S, Q ), (3) 

where Mqq = and Msq = and for consistency (with respect to the heat capacity results), the roots 

of the Eq. m should coincide with two types of the mentioned phase transitions in the heat capacity. It is easy to 
find that due to the structure of the Eq. ©, only in special case Msq = 0 and nonzero Mqq, the divergence points 
of the heat capacity coincide with divergencies of the Weinhold Ricci scalar. In order to obtain consistent results for 
the type one phase transition, the following fine tuning must be hold 

MssMqq - M|q = Ms- (4) 

Regarding various case studies and calculating M and its derivatives, we should note that Eq. is not always 

Af ^ 

satisfied. In addition, it is evident that for the case of Mss = there will be extra divergencies for TZw which 

are not related to any phase transition of the heat capacity. Therefore, the structure of this part of denominator is in 
a way that may present extra divergencies that do not coincide with any type of phase transition points of the heat 
capacity. 
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IV. RUPPEINER METRIC 

The Ruppeiner metric is defined as 0 

d4 = g!!,dY‘^dY\ (5) 

where = —d^S {Y‘^) jdY°'dY^ and = Y°' (M, N'')- It was proved that Weinhold and Ruppeiner metrics are 
related to each other by a Legendre transformation [1^ 

dsji = -MT-^gZdX‘^dX\ ( 6 ) 

Applying this transformation, one can find the following relation for the denominator of the Ricci scalar for this 
case 


denomiUR) = {MssXIqq - T{S, Q)XI^ (S, Q). (7) 

This equation shows that the phase transitions of the heat capacity coincide with the singularities of Ruppeiner 
Ricci scalar only when MJq = 0 and nonzero Mqq. In addition, regarding type one phase transitions, one finds 
due to existence of the T{S,Q) and the fact that T{S,Q) = Ms, the roots of the heat capacity and some of the 
divergence points of Ruppeiner Ricci scalar coincide. On the other hand, there may be some divergence points of 
TZr that do not coincide with phase transition points. These extra phase transition points are originated from the 
zeroes of {MssMqq — M^q). As we mentioned before, in the case of M|q = 0, type two phase transitions of the 
heat capacity {Mss = 0) are covered by divergencies of the Ricci scalar of the Ruppeiner metric. But in this case, we 
encounter with extra divergencies related to the roots of Mqq = 0 which are not related to any phase transition of 

the heat capacity. In addition, in the case of nonzero M^q, the possible real roots of Mqq = lead to the same 
extra divergencies, which were observed in Weinhold metric. 


V. QUEVEDO METRICS 


In order to remove the failures of the Weinhold and Ruppeiner metrics, Quevedo proposed 
metric. The Quevedo metric can be written as [l^ 

a new thermodynamical 

dsl = n {-MssdS^ + MQQdQ^) , 

(8) 

where the (conformal) function Q has one of the following forms 


„ J SMs + QMq, easel 
\ SMs, case II • 

(9) 

In order to obtain the curvature singularity of the Quevedo metric, we calculate the Ricci scalar. Although calcu¬ 
lation of the Ricci scalar is straightforward, analytically calculated results are too large. So, for the sake of brevity 
we do not write the long equations of the Ricci scalar; instead, we can use numerical analysis and some plots to 
investigate the Ricci scalar’s behavior. In addition, since we are looking for the divergence points of the Ricci scalars, 
we can study the denominator of the Ricci scalars with the following explicit forms 

denom{TZQi) = {SMs + QMq)^ MggMQQ, 

(10) 

denom{TlQ2) = S^MgMggMQQ. 

(11) 


Due to Mss being in the denominator of both Quevedo Ricci scalars, the divergencies of the heat capacity and 
Quevedo Ricci scalars coincide. Regarding type one phase transition, although the roots of the heat capacity and 
divergencies of the Quevedo Ricci scalar of case II coincide, for case I this coincidence takes place only for vanishing 
Mq (which is in general a nonzero function). It is worthwhile to mention the fact that there exists an additional 
function Mqq, and its roots provides extra singular points for both cases of Quevedo Ricci scalars. Although Mq and 
Ms are independent from each other, generally, for a nonzero Mq in the case I, it may be possible to set Ms = —^Mq 
for special choices of free parameters. In this situation, one finds another divergence point which may not coincide 
with any phase transition points of the heat capacity. 
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VI. NEW METRIC 

In order to avoid extra singular points in TRS which do not coincide with phase transitions of any type, and also to 
ensure all divergencies of the TRS coincide with phase transition points of the both types, we introduce the following 
new thermodynamical metric 


dslew = {-MssdS'^ + MqQdQ'^) . ( 12 ) 

It is worthwhile to mention that the new thermodynamical metric is defined the same as the Quevedo metric with 
different conformal function. In new metric, we have considered the total mass as thermodynamical potential with 
entropy and electric charge as extensive parameters. In order to hnd the geometrical behavior of new thermodynam¬ 
ical metric, we calculate the Ricci scalar. It is a straightforward calculation to show that the numerator and the 
denominator of TRS for this new metric is, respectively. 


num{TZ) = 6 S'^M|MqqM|_ 5 Mqqqq - QSMIMqqMssMssqq + 25M|qqM|MqqMss 


-f2 


-b4 


SMgMsss - -^Mss {SMss - Ms) 


1 


-MsqMss + MsMssq 


^MqqMsMsqq — 9S^MqqqMsMss 
MqqMsMqqq + '[S'^MgMssQ — S'^MsqMssMsMssq 


SMqqMs {SMss - Ms) Mss - 2 {S'^Mgg + MgMss) Mqq + 2S‘^M‘lQMgg] Mqq, 


(13) 


and 


denom{TV) = S^MgMgg. (14) 

Regarding a typical black hole with nonzero horizon radius, r_|_ (and also nonzero entropy), one can find that, in 
general, the finite mass is an analytic smooth function of its variables. Smooth function is often used technically to 
mean a function that has derivatives of all orders everywhere in its domain. Therefore, one may take into account 
that M has regular derivatives of all orders with respect to the extensive parameters [iq . 

Taking into account the regular numerator of TRS in Eq. (TT^ with the denominator of TRS, Eq. (HI , we ensure 
that all the phase transition points of the type one and two coincide with divergencies of the mentioned TRS and 
there is no extra term that may provide extra divergencies. 

We should note that, in general, the derivatives of all orders of M (such as Mqq, Mgs, Mgq, Mggq and so on) 
are independent from each other. In addition, we should mention that for the case of vanishing Mg or Mgs, the 
numerator of TRS has nonzero value, but denominator of TRS vanishes. In the case of Alg = Mgs = 0, although it 
is clear that both numerator and denominator vanish, the denominator approaches zero faster than the numerator. 
Therefore, one concludes that when Mg and/or Mgs go to zero, TRS diverges. 

In order to elaborate the efficiency of the newly proposed metric and shortcomings of the previously proposed 
metrics, we study two cases of the BI black holes in three and four dimensions. We also discuss linear electrodynamics 
cases and give a comment for neutral solutions. 


VII. BLACK HOLE SOLUTIONS OF EINSTEIN GRAVITY WITH A BI SOURCE 

The d-dimensional action of Einstein-BI gravity in the presence of cosmological constant is given by M 

/ = - / d^xV^ [R-2A + L{F )], (15) 

where R is the Ricci scalar and A refers to the (negative) cosmological constant. Also, L{F) is the Lagrangian of BI 
field as follows 


L{F) = 4/3 




(16) 


where /3 is called the nonlinearity parameter, the Maxwell invariant F = in which F^„ = is the 

electromagnetic field tensor and is the gauge potential. The static BI black hole solutions can be obtained with 
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the following d-dimensional metric 


= -f{r)dt^ + -I- r^dfl 

fir) 


(17) 


where represents the line element of r = constant and t = constant hypersurfaces with volume Vd- 2 - For three 
and four dimensional (3D and 4D) spacetimes, Vd _2 is equal to 2tt and 47 r, respectively, and dOf can be written with 
the following explicit forms 


dV? = 


d6»2, 3D 
dO^ + sin^ Odif'^, AD 

The consistent metric functions /(r) for two cases of three and four dimensional black holes are given by [I9l.l21l 


(18) 


fir) = 


—m — Ar^ -I- 2r^/3^ (1 — T) -|- 


1 - 21 n 


(=AP) 


3D 


_ m _ 


Ar^ 




(19) 


where F = + ■:p 2 {d- 2 )p 2 ; 2 -Fi is hypergeometric function, m and q are integration constants which are related to 

mass parameter and the electric charge of the black holes, resp ectively. The entropy and the electric charge of the 
mentioned B1 black hole solutions were obtained before [13, l2ll - [^ 


5 = 


Vd-2 


r 


d-2 

+ 


4 


( 20 ) 


Vd-2 9 




47r 


( 21 ) 


where denotes the outer (event) horizon of black holes which is the largest real positive root of metric function, 
fir)\r=r+ = 0 . 

Now, we write the quasi-local mass fper unit volume Vd- 2 ) as a function of extensive parameters (S and Q) to 
discuss phase transition. One finds [l^, [2li - l^ 


M 


id-2) 
IGtt 


4(77^Q^-AS^)-|-8S^/3^(1-W)-8ii-^Q^ ln(^(l+'W)) ^ n 


35(7i--A5)+2S^/3^(l-'H)-l-47i-^(3^g 


4D 


( 22 ) 


where ^ = D (^[ 5 ,1] , [|] , — 52^2 ^ and TL = \J4 + Regarding Eq. ([T]), one finds that the heat capacity for 

three and four dimensional B1 black holes can be written as 


{ S''(H-W)[A/3=S=(l-fH)-|-7r''Q''(A-|-2/3''w)l ,, ^ 

2/3iiS(AS^-7r^Q^)(l-|-W)+7i-^Q^A(2+W) 

(23) 

2,32-H^S=*{[E+-l-4(327r=(g-|-4Q25')]W+2Q27i-2-|-6/3^S^} ^ ^ 

~ /3iiW=iS^[E_-l-4Q^7i-^(35+32Qiiy'-|-16Q*g")]W+2Q4.,r4-6/3iiS^(/3iiS^-l-2(5^7rii) 


where = — 35 ' 7 r± 3S‘^{A — 2/3^). It is notable that the prime and double prime are the first and second derivatives 
with respect to Q^, respectively. 

Here, we regard Eq. (1^ and also all obtained TRS in various mentioned methods ( Eqs. ©, (0, ([5|) and (fl^ ) to 
investigate their coincidences. To do so, we plot some figures for three and four dimensional B1 black hole solutions (see 
Figs. wm . We find that for 4-dimensional B1 solutions the Ruppeiner, the Weinhold, both methods of the Quevedo 
and also our new introduced metric have divergencies where the heat capacity diverges (see Fig. [1] and left panel of 
Fig. [3] for more details). In other words, applying the mentioned methods for these solutions leads to coincidence 
between divergencies of TRS and phase transitions of the type two, and therefore, we conclude that for these solutions 
(4-dimensional BI black holes) extra terms in Quevedo metrics, have no contribution to the divergencies of TRS. 

Next, we take into account three dimensional solutions. We find that for all methods of the Weinhold, the Ruppeiner 
and both cases of the Quevedo metrics (Eqs. (I2|), (0 and (I5|)), there are extra divergencies in plotted figures of TRS 


















6 



FIG. 1: BI case: TZ and Cq versus r+ for / = 1, A = —1, j3 = 1, d — 4 and q = 0.1. 

TRS (continuous line) and heat capacity (bold line) for the Weinhold metric (left-up panel), the Ruppeiner metric (right-up 
panel), the Quevedo metric for case I (left-down panel) and the Quevedo metric for case II (right-down panel). 

which come from the contributions of extra term (see Fig. [2|). While for the TRS of our new metric (Eq. (IT^ ) all 
divergencies coincide with the phase transitions of heat capacity (see Fig. [3] (right panel) for more details). 

Another important property of new metric is the behavior of TRS very close to the phase transition points. As Fig. 
[3 confirms, the sign of the TRS before and after the divergence point for these two types of the phase transitions is 
different. If the phase transition is related to the vanishing heat capacity (type one), we see a change of sign for TRS 
before and after of the corresponding singular point. While for the divergence point of the heat capacity (type two), 
TRS has the same sign for the left and right sides of this point. Therefore, this characteristic behavior enables one to 
distinguish these two types of phase transitions from each other. 


A. Linear Case: Maxwell solutions 


In order to study the effect of electric charge and remove the influence of nonlinearity parameter, we investigate 
the Maxwell solutions. We show that these linear solutions elaborate the efficiency of our new metric. In order to 
obtain Maxwell solutions, one can use series expansion of BI solutions for large values of nonlinearity parameter /3. 
Regarding the mentioned results of BI solutions with straightforward series expansion, one can obtain 


—m — Ar^ — In (j) 3D 

4 £) 

r 3 


(24) 


Since the entropy and electric charge of these black holes do not depend on the nonlinearity parameter, one can 
obtain the same results. Therefore, the finite mass (per unit volume Vd- 2 ) of the Einstein-Maxwell solutions can be 
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FIG. 2: BI case: TZ and Cq versus r+ for / = 1, A = —1, j3 = 1, d — 3 and q = 0.1. 

TRS (continuous line) and heat capacity (bold line) for the Weinhold metric (left-up panel), the Ruppeiner metric (right-up 
panel), the Quevedo metric for case I (left-down panel) and the Quevedo metric for case II (right-down panel). 


written as a function of S and Q with the following explicit form 


M = - - - X 

IGtt 






TT 


.2^2 


6U 




Using Eqs. © and (1^51) . one can find the heat capacity may be calculated 


Cq= 


(77=Q"-HAS=)S 

TI-2Q2_AS2 

2S(7rS-AS^-n^Q^) , ^ 

AS2+7I-S-37I-2Q2 


(25) 


(26) 


Using Eq. (1^51) . we are in a position to study the behavior of TRS for different methods that were mentioned in 
this paper. We can use Eqs. @, ®, ® and (US with the heat capacity relations of Einstein-Maxwell solutions, Eq. 
(HH): to plot various figures for studying the geometrical behavior of different thermodynamical spacetime (see Figs. 

H-ED. 

Regarding Figs. |3]-[71 we find that for 3-dimensional Einstein-Maxwell solutions (with considered values for different 
parameters), there is only one root for the heat capacity. In other words, in this case, we see one phase transition of 
type one and there is no divergency for the heat capacity. Whereas in 4-dimensional solutions, there are one phase 
transition of type one and two phase transition of type two. It means that there are three phase transition points in 
four dimensional Einstein-Maxwell black hole solutions. 
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FIG. 3: BI case: TZ and Cq versus r+ for Z = 1, A = —1, 13 = 1, q = 0.1. 

TRS (continuous line) and heat capacity (bold line) for the new metric with d = 4 (left panel) and d — 3 (right panel) 


Regarding Weinhold approach, one finds both mismatching and extra divergencies in three and four dimensions 
(see Figs. |3]and|n|). In other words, in 4-dimensional black holes (Fig. 2]), the root and one of the divergencies of the 
heat capacity are not compatible with any divergency of the Weinhold’s TRS. On the other hand, in 3-dimensional 
solutions (Fig. |6l), this approach completely fails to provide an effective machinery for studying the phase transitions. 
One divergency is observed which does not coincide with any phase transition point and for the root of heat capacity 
TRS is a smooth function. 

As for the Ruppeiner metric, in 4-dimensional Einstein-Maxwell solutions (Fig. 2]), one of the divergencies of the 
heat capacity is not matched with any divergency of the Ruppeiner’s TRS. For 3-dimensional case, there are two 
divergencies for TRS (Fig. 2]). One of them coincides with the root of heat capacity whereas the other one is not 
related to any phase transition point. 

Now we investigate both types of Quevedo metrics and their behaviors. Regarding case I in 4-dimensional solutions 
(Fig. 2] up panels), there is no divergency for TRS in the place of heat capacity vanishing point, whereas divergencies 
of the heat capacity are matched with divergencies of TRS. For case II (Fig. [5] down panels), we obtain effective 
results. In other words, all divergence points and vanishing point of the heat capacity are matched with divergencies 
of the Quevedo’s TRS of case II. But for charged BTZ black holes, regardless of case I or II, an extra divergency is 
seen which is not related to any phase transition point of the heat capacity (Fig. |6]down panels). 

Finally, for both 3 and 4 dimensional Einstein-Maxwell black holes (see Eig. El), new metric is completely successful 
for describing phase transition points of the heat capacity in context of geometrothermodynamics. In other words, 
regardless of the dimensions, the phase transition points of the heat capacity and divergencies of the Ricci scalar of 
the constructed spacetime, coincide and this approach is free of extra divergencies for its TRS. 

We should note that this new metric enjoys the characteristic behavior which was mentioned in Einstein-BI solutions 
and one can distinguish two types of phase transitions from each other. 


B. Neutral (uncharged) solutions 

In this section, we study the behavior of the system in case of uncharged solutions {q = 0). Considering the 
uncharged solutions with an extensive parameter (entropy), one can find that the TRS of mentioned methods of 
geometrothermodynamics vanishes. Therefore, it seems that it is not possible to use these methods for studying the 
behavior of phase transitions. In other words, the geometrothermodynamics methods are valid for systems containing 
two or more extensive parameters. 


VIII. GENERALIZATION 

Now, we are in a position to generalize the new defined metric to the case of more than two extensive parameters. 
Regarding Weinhold, Ruppeiner and Quevedo approaches, one finds that supplementing additional extensive param¬ 
eters increases the complexity of the thermodynamical Ricci scalars and also their denominators. The supplemental 
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FIG. 4: Maxwell case: TZ and Cq versus r+ for Z = 1, A = —1, d = 4 and q = 0.1. 

TRS (continuous line) and heat capacity (bold line) for the Weinhold metric (up panels) and the Ruppeiner metric (down 
panels). 


extensive parameters lead to extra terms which may contribute to the additional number of divergencies of the system 
under study. Considering the total mass of the system as a function of arbitrary number of extensive parameters Xi’s, 
one finds that the denominator of TRS for Weinhold, Ruppeiner and Quevedo approaches will be more complicated. 
Numerical calculations show that denominator of TRS for the mentioned geometrical approaches contain extra terms 
which may contribute to number of the TRS divergencies. For more clarifications, in addition to S and Q, we consider 
angular momentum, J, as extensive parameter. Regarding the total mass of the black holes as a function of S', Q and 
J, we calculate the denominator of TRS for the Weinhold, the Ruppeiner and the Quevedo metrics 


where 


denom{TZ) = < 


{SMs + QMq + JMjf 


Weinhold 
Ruppeiner 
Quevedo Case I 
Quevedo Case II 


i = [Mss {Mqj - MqqMjj) + M^qMjj + MqjMqq - 2MsqMsjMqj\ . 


(27) 


Considering Eq. (EZl), one finds that although vanishing points of heat capacity coincide with related divergence 
points of TRS in the Ruppeiner case (due to the existence of T), for the Weinhold one the coincidences are not 
generally observed. Regarding the divergence points of the heat capacity, which are related to vanishing Mss-, a-nd in 
order to match both divergence points of heat capacity and those of the Weinhold and the Ruppeiner cases, one should 
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FIG. 5: Maxwell case: TZ and Cq versus r+ for Z = 1, A = —1, d = 4 and q = 0.1. 

TRS (continuous line) and heat capacity (bold line) for the Quevedo metric for case I (up panels) and the Quevedo metric for 
case II (down panels). 


set last three terms to zero and the coefficient of Mss should be a nonzero finite expression {Mqj — MqqMjj ^ 0). 
Since these fine tuning conditions are not hold in general, we encounter mismatch between divergence points and/or 
extra divergencies in TRS. 

As for Quevedo metrics, due to existence of M|g, divergencies of both heat capacity and Ricci scalar match with 
each other. As for the phase transition type one and in order to have coincidence between roots of the heat capacity 
and divergence points of the TRS, for the case II Quevedo metric, existence of M| is sufficient to ensure the mentioned 
matching whereas we have the following restriction for the case I 


SMs + QM + JMj = Ms. 


(28) 


In general Eq. (1281) is not hold for arbitrary metric function parameters and therefore, there will be extra divergencies 
for the Ricci scalar which are not matched with any phase transition point. Regarding both of Quevedo metrics, Mqq 
and Mjj, may contribute to additional divergencies of the Ricci scalar which was seen in case of 3-dimensional 
solutions. 

In order to solve these problems we generalize the new introduced thermodynamical metric to the case of n extensive 
parameters (n > 2) with following form 


dSj\[ew — 


SMs 


/tt" O^m Y 
V dx't ) 




(29) 


where Xi Y order to obtain the curvature singularity of the new generalized metric, we should calculate the 

Ricci scalar. Since analytical calculations are too large, for the sake of brevity, we study the denominator of the 
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FIG. 6: Maxwell case: TZ and Cq versus r+ for Z = 1, A = —1, d = 3 and q = 0.1. 

TRS (continuous line) and heat capacity (bold line) for the Weinhold metric (left-up panel), the Ruppeiner metric (right-up 
panel), the Quevedo metric for case I (left-down panel) and the Quevedo metric for case II (right-down panel). 


Ricci scalar. Calculations show that the too long expression of the numerator of TRS (generalization of Eq. (IT^ ) is 
divergence free [l^, while its denominator is 


denom{TV) = S^MgMgg. (30) 

Eq. (1501) confirms that all singular points of TRS coincide with both types of heat capacity phase transition points 
without any extra divergency. 


IX. CONCLUSIONS 

Motivated by a surge of study of geometrical concept of black hole thermodynamics, we introduced a new metric 
regarding the matter. Considering the various approaches toward geometrical study of black hole thermodynamics, we 
first discussed the shortcomings of the mentioned methods. It was believed that the divergencies of TRS indicate the 
thermodynamical phase transitions, and therefore we focused on the roots of denominator of TRS (since we regarded 
smooth function of mass and its derivatives, the numerator of TRS is a regular function and therefore, divergencies of 
TRS is equivalent to the roots of denominator of TRS). Taking into account the Weinhold, Ruppeiner and different 
types of Quevedo metrics, we showed that divergencies of TRS may not coincide with roots and divergencies of the 
corresponding heat capacity. 

In order to avoid this problem, we introduced a new metric. We showed that the denominator of its TRS contains 
terms which are only the product of numerator and denominator of the corresponding heat capacity. In other words, 
all divergencies of TRS in this approach coincide with phase transition points of the heat capacity. Moreover, for 
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FIG. 7: Maxwell case: TZ and Cq versus r+ for Z = 1, A = —1, q = 0.1. 

TRS (continuous line) and heat capacity (bold line) for the new metric with d = 4 (left and middle panels for different scales) 
and d = 3 (right panel) 


more clarifications, we regarded two known examples to show the shortcomings of the other previous metrics and the 
efficiency of the new devised metric. 

Next, we generalized these metrics to contain more than two extensive parameters. As it was seen, in the cases of 
Weinhold, Ruppeiner and different types of Quevedo metrics, (denominator of) TRS contained complicated expressions 
and extra terms that may increase the number of TRS divergencies and shift the place of these divergencies in a way 
that they may not coincide with phase transition points arisen from the heat capacity. We showed that the divergencies 
of TRS related to the generalized new introduced metric are compatible to the phase transition points of heat capacity. 

Another important property of the new metric is the different behavior of TRS before and after its divergence 
points. It was seen that the behavior of TRS for divergence points related to two types of the phase transition is 
different. Therefore, considering this approach also enable us to distinguish these two types of phase transition from 
one another. 

Recently, it was seen and proposed that different constants (such as cosmological constant, BI nonlinearity pa¬ 
rameter, Gauss-Bonnet par ameter, Newton constant and etc.) may vary and have contribution to thermodynamical 
structure of the system [^. In other words, in case of black holes, the total mass of the black hole is a function of 
these parameters as extensive parameters. It will be interesting to reconsider these constants as thermodynamical 
variables and modify TRS of the mentioned geometrothermodynamical methods. Although this modification changes 
TRS of various methods, it does not cause inconsistent results for new generalized metric. 

The approach that we introduced in this paper enable one to map the divergence points of its TRS with phase 
transition points without any concern regarding contribution of other terms and extra divergence points. This method 
may also be employed to study phase transition of other non gravitating systems. 
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